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Summary 


‘An important meteorological input into FSCBG (Teske et al., 1992) is net 
radiation index N, a measure of the influence of solar radiation on atmospheric motion. 
This parameter sets several atmospheric stability variables for the subsequent FSCBG 
calculation (including wind speed power law exponent, azimuthal and elevation standard 
deviations, and mixing layer height); its value impacts the prediction of spray cloud 
growth and expansion. Over the years FSCBG documentation and training session 
materials have focused on several field-observation approaches to establishing the value 
of N. Some rules-of-thumb have been suggested that basically end with setting N = 1 
since the USDA Forest Service recommends spraying in the early morning hours. 


Recently, several members of the FSCBG community have questioned whether N 
could be better quantified, based either on a more analytical approach or by incorporating 
an on-site field measurement of incoming solar radiation. This paper summarizes an 
effort at examining the prediction of N by these approaches. 
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Net Radiation Index 


_The stability character of the atmosphere affects any material released into it. 
Pasquill (1961) first proposed a simple technique whereby dispersion estimates could be 
generated in terms of routine meteorological data. This system was subsequently 
modified by Gifford (1961) into a format widely known today as the Pasquill-Gifford 
curves. These curves enable a researcher to estimate the horizontal and vertical 
dispersion characteristics of a plume release, using readily available meteorological 
information. Cramer (1957) made significant contributions to this technique, while 
Turner (1970) assembled sources of information into a substantial resource workbook for 
estimating atmospheric dispersion. 


Traditionally, atmospheric stability is divided into seven stability classes: 


Stability Stability Stability 
Category Effect Parameter 
A very unstable 1 
B unstable ps 
C slightly unstable 3 
D neutral 4 
E slightly stable 5 
F stable 6 
G very stable 7 


Atmospheric stability is dependent primarily upon net radiation and wind speed. In the 
absence of cloud cover, incoming solar radiation (or insolation) during the day is 
dependent only upon solar altitude (the position of the sun above the horizon), which is a 
function of time of day, time of year and location. Turner (1964) developed a technique 
whereby the solar altitude and actual cloud cover could be used to infer what he called the 
net radiation index. This index could then be correlated with wind speed to determine the 
stability parameter and quantify the stability effect. 


The cookbook procedure works like this (following Turner, 1964, and Fulle, 
1976): 


1) Determine the Insolation Class. From field observations at the time of the spray 


mission, the insolation class number I is determined by knowing the solar altitude 
(where 90 deg is overhead in the sky) using the table: 


Solar Altitude Insolation Insolation Class 
Angle A (deg) Number I 
60<A<90 strong 4 
35<A<60 moderate 3 
15<A <35 slight 2 
0<A<15 weak 1 
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2) Correct for Cloud Cover. The following steps are then taken: 


1. 


If the total cloud cover is 10/10 (overcast) and the cloud ceiling is less than 7,000 


- feet, net radiation index N = 0 (daytime or nighttime). 


input. 


For nighttime (between sunset and sunrise) when the total cloud cover is greater 
eo 4/10 and less than 10/10, N = -1; when the total cloud cover is less than 4/10, 


For daytime when the total cloud cover is less than 5/10, net radiation index 
equals the insolation class number, N = I. 


For daytime when the total cloud cover is greater than 5/10, net radiation index is 
obtained by modifying the insolation class number by: 


a. When the cloud ceiling is less than 7,000 feet, N =I - 2. 


b. When the cloud ceiling is greater than 7,000 feet and less than 16,000 feet, 
N=I-1. 


c. When the total cloud cover is 10/10 and the cloud ceiling is greater than 7,000 
feet, N=I- 1. 


d. Otherwise, N = I. 
e. Plus, N must not be less than 1. 


The resulting value of N is entered at the appropriate place in the FSCBG data 
Table entries are then manipulated by FSCBG to determine the necessary 


meteorological parameters (Dumbauld and Bowers, 1983). The tabular relationship 
determining the stability class, for example, is found from the table (Turner, 1964, and 
Fulle, 1976): 


Wind Speed Net Radiation Index N 
(knots) 4 3 2 1 0 -1 -2 


The curve fits used in FSCBG to extract the necessary atmospheric parameters 


compress the N = -2 value into N = -1 (Bjorklund et al., 1984). This procedure appears 
the most workable, and is strongly recommended. 
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Perr ewaws 


Jim Rafferty (private communication) has provided an equation set to recover the 
solar altitude at any location and time on the earth. A simplification of these equations 
defines the following computational procedure: 


> 


T 


latitude (positive north) 
local time in hours 
0.985648 ( Julian Day - 1 ) 


12.0 + 0.12357 sin d - 0.004289 cos d 
+ 0.153809 sin 2d + 0.060783 cos 2d 


279.9348 + d + 1.914827 sin d - 0.079525 cos d 
+ 0.019938 sin 2d - 0.00162 cos 2d 


arcsin ( 0.39785 sin o ) 
15(T-M) 


solar altitude angle 
arcsin ( sin @ sin D + cos 6 cos D cos h ) 


This information has been collected into a QuickBasic computer program for 
recovering an approximate value for net radiation index. A listing of this program is 
provided in the Appendix of this report. This program assumes that each time zone 
covers 15 degrees of longitude, and that daylight occurs between 6 AM and 6 PM. 


It is anticipated that a more complete equation formulation (developed from the 
Real Time Volume Source Model at U. S. Army Dugway Proving Ground) will be 
included in a later version of FSCBG. 
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Solar Radiation 


_A suggested alternate approach for recovering net radiation index is to make use 
of a direct measurement of incoming solar radiation. The applicability of this approach 
will now be discussed (following Geiger, 1965). 


_ Radiation arrives at the earth's surface from the sun, and reflects by several 
physical processes through the atmosphere. A positive net radiation adds heat to the 
surface of the earth. Radiation is usually measured in calories per square centimeter per 
minute (cal/sq cm min), also called langley per minute (ly/min). Units conversion is: 


1 Y= | cal __ ~ 697.4 W = 291.2 Bo 
min cm2min m2 hr ft2 


___ Net radiation may be denoted by the symbol S. If insolation (incoming solar 
radiation) is greater than outgoing (or terrestrial) radiation, the radiation balance is 
positive; if it is less, the radiation balance is negative. The radiation balance consists of 
two radiation streams of different spectral ranges. The first radiation stream is a short- 
wavelength part available only when the sun shines. Radiation reaching the surface of 
the earth consists of that part of direct incoming solar radiation I not reflected by clouds, 
absorbed by the atmosphere, or scattered diffusely, and that part of the nondirectional sky 
radiation H that represents diffusely scattered radiation that has reached the ground and 
provides "daylight" within the visible spectrum. The value of I + H reaching a 
horizontal surface is called global radiation. Part of this radiation is in turn reflected by 
the earth's surface. This short-wavelength reflected radiation R depends on the nature of 
the ground. Dave Miller (private communication) has provided some insight into the 
measurement of the reflection factor within a canopy, referencing work by Monteith 
(1976). 


The second radiation stream is due to incoming long-wavelength radiation. The 
earth's atmosphere contains water vapor, carbon dioxide, and ozone, all of which absorb 
radiation and re-emit it. The long-wavelength atmospheric radiation G is termed counter- 
radiation since it counteracts the terrestrial radiation loss. It occurs both day and night, 
and in fact is somewhat greater during the day since it depends on temperature. 


In addition, the earth's natural surface cover acts as a black body to emit radiation 
through the soil surface by day and night according to the Stefan-Boltzmann fourth- 


power temperature law oT 4. The radiation balance S is therefore given by the equation: 
S = I1+H+G-oT4-R 


The net radiation S is seen to be a balance of several terms, some (such as G and oT 4) 
potentially an order of magnitude larger than the others. Thus, it seems clear that the best 
field measurement to make would be for net radiation. Unfortunately, no currently 
available instrument is likely to possess the skill needed to produce this measurement (by 
being able to sort out the competing radiation terms). Rather, the direct incoming solar 
radiation I has been found to be measurable, and has thus been correlated with stability in 
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the literature. The stability tables have been reworked for incoming solar radiation, either 
in graphical form (Pasquill and Smith, 1983) or tabular form (Williamson and 
Krenmayer, 1980) as: 


Wind Speed Incoming Solar Radiation (ly/min) 
(knots) 0.0-0.1 0.2-0.3 04 05 £40.6-0.7 0.8-0.9 1.0-1.1 1.2 
0-3 4 2 oe #) ie 1 1 1 
4-5 4 3 4 2 Z 2 1 1 
6 4 4 3 3 3 2 Z Pi 
ti 4 4 3 3 3 3 2 ps 
8 4 4 4 3 a 5 2 Z 
9 4 4 4 3 3 3 a Z 
10-11 4 4 4 4 5 3 C 3 
ihe 4 4 4 4 4 3 3 3 


In this table the incoming solar radiation and wind speed would be used to 
determine the stability parameter, which would then be backfit in the Net Radiation Index 
table. The advantage of this approach is that cloud cover, cloud height, location and time 
of day are not needed; the disadvantage is that the recovery of net radiation index may not 
produce a unique value for N. For example, for incoming solar radiation of 0.2 ly/min 
and a wind speed of 6 knots, a stability parameter of 4 is obtained from the Incoming 
Solar Radiation table (above). A value of 4 in the Net Radiation Index table suggests a 
value of N of 1 or 0. Since FSCBG will accept and use non-integer values of N, the 
appropriate data input entry is probably N = 0.5. 


If this approach proves worthwhile, it could be quantified and included in a later 
version of FSCBG as well. 
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Appendix 
Source Listing of NETRAD.BAS 


CLS 

PRINT “Net Radiation Index Determination" 

PRINT “3-9-3939 - A 

PRINT "Enter the following information as decimal values:" 
PRINT 

INPUT "North Latitude (deg): ", P 

PR = 0.0174533*P 


INPUT "Local Time (0.0 to 24.0 hours): ", T 
DN = 0 

IF T > 6.0 AND T < 18.0 THEN DN = 1 

INPUT “Julian Day (1 to 366): ", J 


D = 0.985648* (J-1) 


DR = 0.0174533*D 
SD = SIN(DR) 
CD = COS (DR) 


STD = SIN(2.0*DR) 

CTD = COS (2.0*DR) 

M = 12.0+0.12357*SD-0.004289*CD+0.153809*STD+0 .060783*CTD 

H = 15.0* (T-M) 

HR = 0.0174533*H 

S = 279.9348+D+1.914827*SD-0.079525*CD+0 .019938*STD-0.00162*CTD 
SRe= 0, 0174533*S 

TEM = 0.39785*SIN(SR) 

DD = ATAN(TEM/SQR(1.0-TEM*TEM) ) 


TEM = SIN(PR) *SIN (DD) +COS (PR) *COS (DD) *COS (HR) 

A = 57.29578*ATAN (TEM/SOR(1.0-TEM*TEM) ) 

I= 4 

IF A < 60.0 THEN I = 3 

rhe AwaS 5. 0) THEN L =) 2 

DreAeNel OO THEN T= 1 

INPUT "Cloud Cover Fraction (0.0 to 1.0): ", CC 


CH = 0.0 
IF CC > 0.5 THEN INPUT "Cloud Ceiling in 1000s of feet (1 to 30): ", CH 
IF>cC > 0.95 AND CH < 7.0 THEN 
N = 0 
ELSE 
IF DN = 0 THEN 
IF CC < 0.4 THEN 


Nis 2 
ELSE 
N = -l 
END IF 
ELSE 
Nz=tTI 


LeaCCe>mve. so LHEN 
IF CH < 16.0 THEN N = I-1 
IF CH < 7.0 THEN N = I-2 
IF cc > 0.95 AND CH > 7.0 THEN N = I-1 
IF N < 1 THEN N = 1 
END IF 
END IF 
END IF 
PRINT 
PRINT "Net Radiation Index: "; N 
END 
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